Spin filtering through excited states in double quantum dot pumps 
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Recently it has been shown that ac-driven double quantum dots can act as spin pumps and 
spin filters. By calculating the current through the system for each spin polarization, by means of 
the time evolution of the reduced density matrix in the sequential tunneling regime (Born-Markov 
approximation), we demonstrate that the spin polarization of the current can be controlled by 
tuning the parameters (amplitude and frequency) of the ac field. Importantly, the pumped current 
as a function of the applied frequency presents a series of peaks which are uniquely associated with 
a definite spin polarization. We discuss how excited states participating in the current allow the 
system to behave as a bipolar spin filter by tuning the ac frequency and intensity. We also discuss 
' spin relaxation and decoherence effects in the pumped current and show that measuring the width 

(~ i ' of the current vs frequency peaks allows to determine the spin decoherence time T2 . 
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I. INTRODUCTION 



\ Production, control and detection of spin polarized current through nanostructures (spintronics)^ has become an 
area of intense activity in the past few years. This is mostly due to the long coherence times of spins, as compared to 
charge, giving rise to possible applications in quantum information processing. Proposals for generating spin-polarized 
. currents include spin injection by using ferromagnetic metals^ or magnetic semiconductors^. Alternatively, one may 
S ' use quantum dots (QDs) as spin filters or spin pumps^j^^Sii^. In a semiconductor QD, the number of electrons can 
be controlled via gate voltages, making it an ideal system to study transport of individual electrons and probe the 
structure and properties of the discrete energy level spectrum. Indeed, the spin of an isolated electron in a QD has 
been proposed as a quantum bit for transport of quantum information. For QD spin filters, dc transport through 
few electron states is used to obtain spin-polarized currents which are almost 100% spin-polarized as demonstrated 
1^ [ experimentally by Hanson et af following the proposal of Recher et a? . The basic principle of spin pumps is closely 
related to that of charge pumps. In a charge pump a dc current is generated by combining ac driving with either 
' absence of inversion symmetry in the device, or lack of time-reversal symmetry in the ac signal. The range of possible 
\^ I pumps includes turnstiles^, adiabatic pumps or non-adiabatic pumps based on photon-assisted tunneling (PAT)iS*ii. 
' In the last few years, the application of ac electric fields in quantum dots has shown to be very accurate both to control 
, and to modify their transport and electronic propertiesifi. For instance. Sun et alJ^ have proposed a spin cell based 
on a double quantum dot (DQD) driven by microwave radiation in the presence of an external non-uniform magnetic 
field. For strongly correlated electrons in quantum dots in the Kondo regime, the application of an ac potential has 
been shown to induce decoherence and destroy the Kondo peal*i2*i^. An ac potential modifies as well the electron 
G ' dynamics within single and double quantum dots and it has been shown that under certain conditions it induces 
1 charge localization and destroys the tunneling, allowing, by tuning the parameters of the ac field, to control the 
time evolution of the state occupation, and therefore the entanglement character of the electronic wave function^°-^^. 
O , Regarding detection of single spin states in semiconductor quantum dots, this has been achieve d^^d^ using quantum 
O ' point contacts and spin to charge conversion. These studies allow determination of spin relaxation times on the order 
^ , of milliseconds. 

Recently, we proposed a new scheme for realizing both spin filtering and spin pumping in an ac-driven DQD (PAT 
regime), in the Coulomb blockade transport regime, by tuning the parameters of the ac fieldi^. We considered an 
^ , asymmetric DQD system with one energy state in the left dot and two orbital states in the right dot, allowing 
each QD to contain up to two electrons. An in-plane magnetic field breaks the spin degeneracy producing Zeeman 
splitting. Thus, for the right dot, we can have the following singlet and triplet states, in order of increasing energy: 
1^0) = ^ (lit) - I Ti)), \T+) = I rn, \To) = ^ (I it*) + I Ti*)), |T_) = I ii*) and \S^) = ^ (| it*) - | Ti*)), 
where the electrons in the upper level are marked with an asterisk (*). In our previous work^*, we considered a 
basis up to the lowest energy triplet state |T+), and we discussed the conditions for obtaining spin polarized current 
from unpolarized leads. Here, we complement these results by adding the remaining excited triplet states which 
also participate in the pumping process. We find that at certain frequency of the ac field and in the presence of an 
uniform magnetic field the double quantum dot acts as a pump of fully spin polarized electrons. Increasing the ac 
frequency, we show that an inhomogeneous magnetic field is necessary in order to obtain fully spin polarized current 
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in the reverse direction. Then the system, in such a configuration, acts as a bipolar spin filtering just by tuning the 
ac frequency. 

The periodic variation of the gate potentials in our system, consisting on an asymmetric DQD, allows for a net dc 
current through the device even with no dc voltage applicd^^-^*': if the system is driven at frequencies (or subharmonics) 
corresponding to the energy difference between two time-independent eigenstates related by the intradot tunneling of 
an electron, this electron become completely delocalized^^'^^. If the left reservoir (chemical potential /il) can donate 
electrons to the left dot (at a rate r^) and the right reservoir (chemical potential /i^) can accept electrons from the right 
dot (at a rate r^^) the system will then pump electrons from left to right, even when there is no dc bias applied, namely 
Mi — M-R- Starting from this pumping principle our device has two basic characteristics: i) If the process involves two- 
particle states, the pumped current can be completely spin-polarized even if the contact leads are not spin polarized 
and ii) the pumping can occur cither through triplet (Fig. la) or singlet (Fig. lb) states depending on the applied 
frequency, such that the degree of spin polarization can be tuned by means of the ac field. For example, if one drives 
the system (initially prepared in a state with n = Ul + tir = 3 electrons: \L =['[, R =T)) at a frequency corresponding 
to the energy difference between the singlets in both dots, namely oJSol ~ Esg.R — Eso,l (we consider h = e = 1), the 
electron with spin | becomes delocalized in the DQD system. If now the chemical potential for taking | (t) electrons 
out of the right dot is above (below) ^r, a spin-polarized current is generated. The above conditions for the chemical 
potentials can be achieved by breaking the spin-degeneracy through a Zeeman term in each dot — \g\iiBBa 
(a — L,R), where Ba is the external magnetic field at dot a, (which is applied parallel to the sample in order to 
minimize orbital effects), g is the effective g- factor and fiB the Bohr magneton. Then, for example, a fully spin-down 

polarized pump is reahzed at the frequency utsal through the sequence: (it, T 0) ^ (t, it 0) =P (t, T 0) ^4 (it, T 0) 

or (it,T 0) (T,iT 0) =^ (it, it 0) ^ (iT,T 0) which involves singlet states. Increasing the frequency of the ac 
voltage, the states (J,t, T 0) and {l,T^) arc brought into resonance at wt+t ~ Et+,r — Es^.l and spin up polarized 
current flows to the collector from the right quantum dot. However, it can be seen that the state (|t, T 0) is also in 
resonance with (t,To) at the frequency ujtqI ~ Eto,r — Es„,l — '^t^^ + ^r — A^, which coincides with wt+t for an 
homogeneous magnetic field, A^; = A^, and this process contributes to spin down current. Consequently, the current 
at this frequency is partially spin polarized. This degeneracy, present in the three electron configuration, is removed 
if the Zeeman splitting is different in the left and right dot. This could experimentally be achieved either by applying 
a different magnetic field to each dot or considering quantum dots with different g-factors. In that case, wt+^ ^TqI 
and a pure spin up current is obtained by applying an ac potential with frequency ujt+]- 

In practical implementations of such devices, the effects of spin relaxation and decoherence, characterized respec- 
tively by times Ti and T2, need to be addressed. The most important source of decoherence in quantum dots is the 
hyperfine interaction between electron and nuclear spins. Recently, several studies have been reportedSiSiiS^ measur- 
ing these characteristic times and also T2*, the spin dephasing time for an ensemble of nuclear spins. Decoherence 
effects are included in our model in a phenomenological way and we give an analytical treatment to explain the 
numerical results. Interestingly, we find that the decoherence time T2 can be obtained from the width of the current 
vs frequency peaks. 

This paper is organized as follows: in section II we introduce our model and write the equations for the time 
evolution of the density matrix. In section III we describe the results obtained for homogeneous and inhomogeneous 
magnetic field. In section IV we study spin relaxation effects both numerically and analytically, and discuss the 
possibility of obtaining the decoherence time T2 from the width and height of the current vs frequency peaks. We 
finalize with conclusions and outlook in section V. 



II. THEORETICAL MODEL 

Our system consists of an asymmetric DQD connected to two Fermi-liquid leads which are in equilibrium with 
reservoirs kept at the chemical potentials fj,a, a = L,R. Using a standard tunneling Hamiltonian approach, we 
write for the full Hamiltonian Hi + Hdqd + Hz + Ht, where Hi — J2a a ^^^c]. ^Ck^a describes the leads and 
T^DQD = 'Hqjj + ^QZ5 + "^l^^r describes the DQD. Hqjj,Hqjj describe each dot including the charging energies 
of the dot electrons. The presence of an external magnetic field is taken into account through the term Hz = 
i ^ jj., Actdjj, g.((T2)(jcr'rfa,cr' , where Aq = g^sBa is the Zeeman splitting of the energy levels of each QD in the 
presence of an external magnetic field Ba — (0,0, Bq,). The model has been described in detail elsewhere^* but we 
include the description here for the sake of clarity. It is assumed that only one orbital in the left dot participates in 
the spin-polarized pumping process whereas two orbitals in the right dot (energy separation Ae) are considered. The 
isolated left dot is thus modelled as a one-level Anderson impurity: Hqjj = E1d\^dLa + ULnL'\nLi, whereas the 
isolated right dot is modelled as: H^j^ = J^ia ^m^RiadRia + UR{J2i nRi^UR^i + Y.a,a' "flOo-nflio-') + -^SoSi including 
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the exchange interaction term. The index z = 0, 1 denotes the two levels. In practice, we take e]^ — Ej^^ = 0, 

E^ = Al and E^^q — Ar, with charging energies Ur > Ul- Si = {l/'^)J2aa' d?Ricr'^<y<y'dRia' are the spin operators 
of the two levels. As a consequence of Hund's rule, the intra-dot exchange, J, is ferromagnetic (J < 0) such that 
the energy of the singlet \Si) = (l/\/2)((i^g|d^^ ~ ^AOi'^mT)!'^) i^SuR. = Ur + Ar + Ae — ^) is higher than 
the energy of the triplets \T+) = 4ot4it|0>> (^t+m = Ur + As + {), |To) = (l/x/2)(4o^4i^ + 4oi4iT)|0) 
{Eto,R = Ur + Ar + Ae+ {), and |T_) = 4oi4ulO> (Et^i-.R ^ Ur + 2Ar + Ae + |). As can be seen, due 
to the Zeeman splitting, Et_ > Etq > Et^. Finally, we consider the case where Ae > Ar — J/4 such that the 
triplet |r+) is higher in energy than the singlet |S'o) = (l/V2)(4oT^flOi " d^Roi'^m^l'^) (^SoX(fl) = Ul{r) + Aur,)). 
T^L-^R — X^i (J ^i-R(c^Lcr^-Ri(T + ^-C-) describes interdot tunneling. The tunneling between leads and each QD is described 

by the perturbation Ht — J2kL a- ^L{clj^cr'^La+'i^-c ) + J2i Rr a ^i?(Cfe^o-'^-Rjo'+h.c.). In addition we consider an external 
ac field acting on the dots, such that the single particle energy levels become time dependent: 



Vac 



where a =tii: ^Ac and uj are the amplitude and frequency, respectively, of the applied field. We have considered a 
basis of 40 states in the particle number representation which are obtained under these conditions, considering up to 
two electrons in each QD. To study the dynamics of a system connected to reservoirs one can consider the reduced 
density matrix (RDM) operator, p = trRX, where one traces all the reservoir states in the complete density operator of 
the system, x- The evolution of the system will be given by the Liouville equation: '^^^p — —i[H{t), p{t)]. Assuming 
the Markov approximation^^, we obtain the master equation, written as: 

P{t)m'm — —ii^m'm{t)pm'm{t) — i[HL^R, p{t)]m'm 




^ mkPkk kmPmm v^J I ^m'm 

{t){l-5^,ra) (2) 



where ^^/^(f) ~ E„i'{t) — Em{t) is the energy difference between the states |m) and |m') of the isolated DQD, T^k 
are the transition rates for electrons tunneling through the leads, from state j/c) to state |77i) and ^rn'm describes the 
decoherence of the DQD states due to the interaction with the reservoir. This decoherence rate is related with the 

transition rates by: '^'jra'm = \ (Sfc/m' ^fem' + Y^k^m^kn^j- 

Neglecting the influence of the ac field on tunneling processes through the leads, one can calculate these rates by 
the Fermi Golden Rule approximation: 



Tmn = y^,r;(/(a;mn " Ml )<^Af™ , W„ + l + (1 " /(Wnm " W ))'5Af„ ,Af„ - 1 ) , (3) 



where F/ = 27rI?;|V;p; I = L,R are the tunneling rates for each lead. It is assumed that the density of states in both 
leads T^L R and the tunneling couplings Vl.r are energy- independent. Nk is the number of electrons in the system 
when it is in state |fc). 

We include spin relaxation and decoherence phenomenologically in the corresponding elements of the equation for 
the RDM. Relaxation processes are described by the spin relaxation time Ti = {W^i + VFj-f)"^, where Wn and Wi^ 
are spin-flip relaxation rates fulflUing a detailed balance: W|f = exp{—Az/kBT)W^i, where ks is the Boltzman 
constant and T the temperature. A lower bound for the spin relaxation time Ti on the order of /is at i? « — 2T was 
obtained by Fujisawa et ali^I.. Recently a value of Ti — 2.58 ms with a field B = 0.02 T was measured^^ for a single 
QD using a tunnel-rate-selective readout method. In the following, we focus on zero temperature results such that 
Wii = and thus Ti = W^^]^- is the spin decoherence time, i.e., the time over which a superposition of opposite 
spin states of a single electron remains coherent. Recently, Loss et ali^ obtained that T2 — 2Ti for spin decoherence 
induced by spin-orbit interaction. This time can be affected by spin relaxation and by spin dephasing time T2* ^ i.e., 
the spin decoherence time for an ensemble of spins. For processes involving hyperfine interaction between electron 
and nuclear spins, Petta et al.^'^ have obtained T2* w 10 ns from singlet-triplet spin relaxation studies in a DQD. Here 
we consider two cases: T2*^Q.lTi and T2* = O.OOlTi. 

In practice, we integrate numerically the dynamics of the RDM in the chosen basis. In particular, all the results 
shown in the next paragraphs are obtained by letting the system evolve from the initial state | [] ^ ]) until a stationary 
state is reached. The dynamical behavior of the system is governed by rates which depend on the electrochemical 
potentials of the corresponding transitions. 
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We do not include second order processes as co-tunneling. The current from the right dot to the right contact is 
given by: lL~*R{i) = ^R^s Pssit), with a similar expression for Ir^l- Here, states |s) are such that the right dot is 
occupied. 

III. RESULTS 

We analyze the current through an ac-driven double quantum dot, with Vac and oj the amplitude and frequency 
of the ac-field, weakly coupled to the external leads, in the presence of a magnetic field which induces a spin splitting 
in the discrete states of each dot. We assume the ground state in each dot to be the one with spin-up. The model 
that we propose is of general application, i.e., it is valid for both small or large inter-dot coupling. In particular we 
show results for a configuration where the coupling of the dots with the leads is weak and symmetric: T^ ji — 0.001, 
the hopping between the dots \s Ilr = 0.005 and the charging energy for the left and right dots are Ul — 1-0 and 
Ur = 1.3 respectively. All energy units are in meV. The Zeeman splitting produced by an external homogeneous 
magnetic field is: = A^; = 0.026, corresponding to a magnetic field B « IT; the exchange constant for the 
right dot is J = —0.2 and the chemical potentials in the left and right leads are — I^r — 1.31, respectively. We 
have considered two levels in the right dot with energy separation Ae = 0.45 and we neglect inter-dot Coulomb and 
exchange interactions. 

In order to configure the system in a way that electrons can tunnel from the right dot to the right lead only through 
the transitions: | ti)i?. ^ I T)i?. we choose the energy parameters so that they satisfy: Ur < ixr. < Ur + Ar. The 
energy cost of introducing a second electron with either spin-up or spin-down polarization in the left dot has to be 
smaller than the chemical potential of the left lead. This is always satisfied if /i^ > Ul + A^. We consider the system 
to be in the pumping configuration /i^ = fiR throughout. 

Then, if the DQD is initially in the state | tii T)j no current will flow through the system unless the ac frequency 
fits the energy difference between the different states for the left and right dots. For instance: | ti, T) "^^^ I Tj Ti) (at 
Lo = ujsol = C/i? - C/l + Afl - Ai), I ti, T> ^ I T, To) (at = ujt„i = Ur -Ul + Ar - A^ + Ae + J/4), | Ti, T) ^ I T, ^i) 
(at uj = ujs.i = Ur-Ul + Ar - Al + Ae^ 3J/4) or | ti, T) ^ I i,T+) (at w = wt+t ^Ur-Ul + Ae + J/4). 
The suffix i (i) is used to remark that the inter-dot tunneling electron has spin- up (down) polarization. At these 
frequencies, one electron becomes delocalized undergoing Rabi oscillations with frequencjsifi 

Qr = 2tLRM — ), (4) 

where Ji, is the vth-ovdev Bessel function of the first kind and v is the number of photons required for bringing the 
two states into resonance. It is important to note that the frequencies for transitions involving spin down depend on 
the difference Ar — Al while the one involving spin up does not. As we will show below, this is the main reason 
for requiring an inhomogeneous magnetic field to obtain spin- up polarized current. One should remark also that the 
resonance achieved between | tiii) I T,^^-) occurs at a photon frequency: ujt^i = ^Tol- Similarly, the resonance 
I Ti, i) ^ I i,T()) occurs at wtoT = ^T+p 
We have analyzed both homogeneous (A^ = A^) and inhomogeneous (A^ ^ Al) magnetic field cases. 

l) Ar = Al. 

As has been discussed above, by tuning the ac frequency a series of spin-polarized current peaks are obtained 
which can be identified and related to different inter-dot resonant processes. If the Zeeman splitting is the same in 
both dots, inter-dot tunneling processes involving electrons with different spins, as for instance: | Ti,T) I T,2o) 
(at to — lutqi) and | Ti, T) ^ I i,7+) (at lu = ti^T+t) occur at the same frequency, i.e., ujtoI — ^t+t ~ ^■^'^ and 
we find two overlapping peaks (Fig^a)). The current is created through several processes. The dominant ones are 
I Ti,T) ^ I T,?o) {| T,T)or| U:To)} -> I Ti,T), which contributes to spin-down current, and | Ti, T) ^ I l,T+) 
{| i,T)or| Ti,T+)} I Ti,T), which contributes to spin-up current. So, in this case, the current is partially spin-up 
polarized. 

The frequency uj — ujsqI ~ 0.3 (Fig. ^b)), corresponds to the one-photon resonance | Ti, T) I T, Ti), which is re- 
sponsible of the large spin-down current peak through the sequence: | Ti, T) ^ I T, Ti) ^ {I T, T)or| Ti, Ti)} ^ I Ti, T)- 
These are the only processes that take place at this frequency, so no spin-up current is expected. Two peaks appear 
in the vicinity of i-^SoIj one at w = ujsii/<i « 0.283 (which corresponds to the three-photon satellite of resonance cosii 
(not shown)) and two overlapping peaks at w « 0.325 (corresponding to the two-photon satellites at ujtoi/'^ ~ 
see Fig^a)). The positions of these peaks are completely independent of each other and are determined by the 
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FIG. 1: Pumped current as a function of the ac frequency (Al = Aii), for the resonances (a) wtqI = i-OT_f_i ~ 0.65 and (b) 
uJSol ~ 0.3. At this frequency, the current is fully spin down polarized. The smaller peaks in (b) are multi-photon satellites 
of other processes: ai lo ~ 0.283, the three-photon process corresponding to the resonance between the singlet So in the left 
dot and the singlet Si in the right dot occurs. The two overlapping peaks at lj « 0.325 correspond to two-photon satellites 
of the resonance in (a). Parameters: Ilr = 0.005, F = 0.001, Ul = 1.0, Ur = 1.3, Ai = = 0.026, Ae = 0.4, J = -0.2, 

fJ-L = fJ-R ^ 1.31, Vac = LOT+T- 



energetics of the system. Thus, we obtain fully spin-down polarized current at lu ~ (^Sol 
filter for spin-down electrons. 



i.e., our device acts as a 



a) Arj^Al. 



In order to get a fully spin-up polarized current peak, we consider different Zeeman splittings between both dots. 
This introduces a separation Au « (where n is the number of photons involved in the resonant transition) 

between peaks with different spin polarization, as can be seen for example by comparing figures ^a) and[2Ia). 
Similarly, the resonances u;Toii= '^T_i)are shifted by an amount A^ — with respect to the resonance at cjt+t(~ 
uJTot) which is independent of Zeeman splitting (see Figl^Ja)). So we obtain fully spin-up polarized current at 
t^T+T ~ 0-65 and fully spin-down polarized current at ujtqI ~ 0.676. In Figl^Jb), the resonance uj = lusqI ~ 0.326 
is not well resolved because there are overlapping satellite peaks (at u — uiTni/^ ~ 0.338, u = wt+t/^ ~ 0.325, 
UJ = LOSii/'^ ~ 0.292). In concrete, we find different processes that contribute to opposite spin polarization currents 
and depend on the absorption of a different number of photons (therefore, their Rabi frequencies are renormalized 
with Bessel functions of different orders, Eq. |0J). It has been shown that for certain ac parameters and sample 
configurations^^the height of the current peaks depend on flj^, which is a non linear function on the ac intensity (see 
Eq. Q). Thus, we can manipulate the spin polarization by tuning the intensity of the ac field. In Fig. |31we show the 
current and current polarization as a function of ac-field intensity at lo=ujsoI- Here, the spin-up contribution comes 
from a two-photon resonance (uj = wt+t/^ ~ 0.325) and, thus, is expected to vanish when J2{c^) = due to the 
dynamical localization phenomenal^. Furthermore, for ac intensities such that Ji(a) = 0, the one-photon resonance 
Losoi quenched and we obtain at this frequency fully spin-up polarized current. Thus, if we tune the ac intensity 
to the value where the first (second)-order Bessel function vanishes, we obtain fully spin-up(down) current. Fig|3|[a) 
shows that at low ac intensities, the contribution of multi-photon processes is small and the tusol ~ 0.326 resonance 
corresponding to practically fully spin-down current is clearly resolved. 



IV. SPIN RELAXATION EFFECTS 



It is important to note also that, contrary to the case for spin-down pumping, the pumping of spin up electrons 
leaves the double dot in the excited state | i, t)- This makes the spin- up current sensitive to spin relaxation processes. 
If the spin | decays before the next electron enters into the left dot, a spin-down current appears through the cycle 
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FIG. 2: Pumped current as a function of the ac frequency (Al 7^ ^r), for the resonances (a) o^TqT = ^t^i ~ 0.65 and 
ujtqi ~ i^T_j^ ~ 0.676 and (b) i^Sgi ~ 0.326. The current at the resonance uj — ujsqI is partially spin down polarized due to 
the overlapping satellite peaks (at lo — ujToi/2 ~ 0.338, lo = lot^i /2 ~ 0.325, uj = ujSii/3 ~ 0.292). Parameters: tLR = 0.005, 
r = 0.001, Ul = 1.0, Ur = 1.3, Al = 0.026, Ar = 2Al, Ae = 0.4, J = -0.2, Vac = c^t+t- 




FIG. 3: (a) Current and (b)current polarization (defined as P = (J| — /|)/(/| + I^)) dependence on the ac field intensity, for 
ujsgi = 0.326. There is fully spin-up polarization for intensities such that Ji{a) ~ and fully spin-down-polarization when 
J2(a) ~ 0. The sample parameters are the same as in fig. 2. 

(it; T) ^ (ii =^ ill T) (T: T) =^ (Ti: T) and the pumping cycle is no longer 100% spin-up polarized leading to 
a reduction of the spin-up current at this frequency. 

We include a finite spin-flip relaxation probability, = l/Ti = 5 x lO^^meV and T2*=0.1Ti and calculate the 
current at ojt+t = '-'•65 for different values of the coupling with the contacts, F, at Vac = 0.14 meV (Fig. ^J. The full 
widths (FWHM) of the current peaks are plotted as a function of F in Fig. I3a) for weak (circles) and strong (squares) 
AC- field intensity. Vac- In order to minimize nonlinear effects, in Fig. EJb) we investigate the low intensity regime 
{Vac = 0.14) where we expect a FWHM dominated by decohercnce, for two different values of the spin dephasing 
time: T2* = G.lTi and T2* ~ .OOlTi. In every case, we find that for large F, the behavior of the FWHM is linear 
with a slope which approaches 2, i.e., FWHM ~ 2F, which can be directly related to the decohercnce time T2 as we 
show below. Thus, experiments along these lines would complement the information about decoherence extracted 
from other setups^S. As an illustration we present below an analytical treatment which allows us to relate the current 
peak widths with the spin decoherence time. 
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FIG. 4: Pumped current near resonance lo = 0.65 for different symmetric couplings to the leads. = A^/0.3, Wn = 5 x 10 ®, 
Vac = 0.14 and T2* = O.lTi. 
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FIG. 5: Full width at half maximum of the total current in Fig. ^ as a function of coupling to the leads, F, at frequency 
ui = 0.65 (a) for strong (squares) and weak (circles) field intensities and T2* = O.lTi and (b) for different spin dephasing times 
= .ITi (circles) and T2'' = .OOlTi (squares), for the weak field intensity case. In both cases, the numerical results are 
compared to the analytic prediction (dashed lines) given by I12I I . (a) In the weak field case, Vac ~ 0.14, Qr « 0.001 (circles), 
we see that for F > ^lR, the curve follows a linear behavior: FWHM 2^ = 2(W||/2 + I/T2* + F). For the high field intensity 
case yAC=0.7, fijj ~ 0.004 and the same behavior is expected at larger F. In (b), for T2* = .OOlTi (squares), we see that the 
width of the peak is larger than the expected asymptotic behavior for large F. This is due to overlapping with the peak at 
i^Toi- As discussed in the text, the extrapolation of the asymptotic curves at F = gives the value of 2(l/2Ti + 1/T2*). Thus, 
for the case where T2* <C Ti, would allow to estimate the value of the spin dephasing time. In both graphics, = A^/0.3, 
Wu = 1/Ti = 5 X 10"'' . 



A. Analytical treatment 



In the following we present an analytical treatment in the stationary regime of the case discussed above where the 
influence of spin-flip processes on the spin up current peak coming from the ac-induced resonance between |1) = | Tii T) 
and |2) = | |, 11*) was numerically obtained. 

As discussed in the previous section, at w = these states are brought into resonance and the current, for 
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An « Al/0.3, is fully spin-up polarized. Intermediate states are |3) = | i,T), |4) = | Ti,TT*), |5) = | T,TT*) and 
|6) = I TiT)- Fro m^?- ™ it is known that the dynamics of the system at large time scales is obtained by a time- 
dependent basis transformation on the density matrix (rotating wave approximation, RWA) such that e{t) — > eo — nui 
and t^n —^ t^n = (— 1)" J„(VAc/^)^Li?,- The equations of motion for the corresponding reduced density matrix 
elements in the RWA are: 



Pi = -2tLR^P2.1 + TlAs + r_R/34 -I- TlPg 
P2 = 2ii,fl9p2,i - (r_R, + Tl + VK|x)p2 



(5) 



while for the intermediate states, 



P3 = ^FI.P2 - (Fl + W^i)p3 

Pi = ^LP2 - ^RPi + ^LP5 
P5 = WiiP2 - {Tl + ^r)P5 
P6 = -^LP6 + ^RP5 + W^iP3 

The equation for the off-diagonal density matrix element is 

P2.1 = [«(eo - nuj) + 7]P24 + iiLR{p2 - 
where 7 = I/T2 is the decoherence rate: 



Pi) 



1 



7 ^ 



(rz 



W^ 



TiJ 



1 



(6) 
(7) 

(8) 



This, together with the condition of conservation of probability SiP; = 1 gives for the total current (at this frequency), 
in the stationary regime, an expression which we can writ*^" as: 



+ (eo - nuY 



where /q = 2^t\J^/W'^ is the current maximum and W is the half width at half maximum: 



n 



n 



+ 7^ 



(9) 



(10) 



Here, F = (F^ -|- Y jV n. For the symmetric case, F^ = Fj^ = F (f = 4) we can rewrite this in terms of the Rabi 
frequency fi/j = Ulr as : 



1-K 



n 



+ 7^ 



(11) 



In our calculations, we have taken 
(fTT|l simplifies to 



lY^W^l \ 8F 
5 X 10^^ ^ ^^7?, F throughout the range of values considered. Then, Eq, 



In the limit W^i,l/T^ < F, then 7 w F, and we get VF^ ^ ^2^ 
result£i9-20_ 

From H12() . we obtain the following asymptotic behaviors: 

1. F <C ilfl; (strong inter-dot tunneling)=> VF^ w /T 

2. F ^ VIr (weak inter-dot tunneling) TF^ w 7 

3. F « ^W^^ 7fl(^^fl + 1r), where 7^ = ^r + I/T2*. 



(12) 



F^ in agreement with previous analytical 
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In Fig. 13 a) the full width at half maximum (FWHM) of the spin-up current peak is represented as a function of F 
for the case T2* — O.lTi together with the analytical curve 1)11(1 . For the weak field case Vac — 0-14 (full circles) with 
riij « 0.001, we see that the predictions of the theory are indeed fulfilled. In particular, in the range F 3> ilfl, the 
FWHM as a function of F is a straight line with slope 2 as expected (FWHM fa 27 = 2/T2* + l/Ti + 2F). This means 
that a direct measure of the decoherence time for the isolated system Tj^": l/r|'*°=l/r2* + l/2Ti can be obtained 
from this linear behavior. From Fig. El^a) it can be verified that the cases (i) {W"^ « a + 6/(rT'2'*°)) and (iii) arc also 
reproduced. Besides, it is interesting to mention that from these cases one can get information on the Rabi frequency. 
The same analysis also holds for the strong field case Vac = 0.7 (full squares), where w 0.004. 

Recently, it has been measured that the spin dephasing time T2* induced by hyperfine interaction is tens of 
nanoseconds^^iS^. This, together with experimental values for Ti as long as miliseconds^^ in GaAs quantum dots, 
such that 1/Ti <C l/?2*, would allow to estimate the spin dephasing time T2* directly from the intersection of the 
large F asymptote with the vertical axis (Fig. |SJb), dotted line). 

We note, from Fig. I3b), that the numerical results differ from the expected analytical curves for the case T2* — 
lO^'^Ti for large F. This is because the parameters that contribute to FWHM (F or 1/12*) are large enough to make 
the current peak overlap with its neighboring spin-down current peak at wtqI- Then, it loses its lorenztian shape, 
mixes its spin polarization and gets wider than what it is analytically expected. This problem should hold in the 
experimental measurements unless the energy difference between both peaks were large enough so they do not overlap 
and they can be fitted to a lorenztian curve. This is the case when — ^ 27. 



V. SUMMARY 



We have proposed and analyzed a new scheme of realizing both spin filtering and spin pumping by using ac-driven 
double quantum dots coupled to unpolarized leads. Our results demonstrate that the spin polarization of the 
current can be manipulated (including fully reversing) by tuning the parameters of the ac field. For homogeneous 
magnetic field, — A/j, we obtain spin down polarized current involving singlet states in both dots. In order to 
obtain spin up polarized current, an inhomogeneous magnetic field is required to break the degeneracy in transitions 
involving triplet states in the right dot. Our results also show, both analytically and numerically, that the width in 
frequency of the spin-up pumped current gives information about the spin decoherence time T2 and also about the 
spin dephasing time T2* of the isolated double quantum dot system. Experiments along these lines would allow to get 
information, from transport measurements, on the different mechanisms producing spin decoherence in quantum dots. 

We thank J. A. Maytorena for useful discussions. Work supported by Programa de Cooperacion Bilateral 
CSIC-CONACYT, by the EU grant HPRN-CT-2000-00144 and by the Ministerio de Ciencia y Tecnologfa of Spain 
through grant MAT2002-02465. 



^ S. A. Wolf, D. D. Awschalom, R. A. Buhrman, J. M. Daughton, S. von Molnar, M. L. Roukes, A. Y. Chtchelkanova, and D. 

M. Treger, Science 294, 1488 (2001). 
^ M. Johnson and R. H. Silsbee, Phys. Rev. Lett. 55, 1790 (1985); F. J. Jedema, A. T. Filip and B. J. van Wees, Nature(London) 

410, 345 (2001). 

^ R. Fiederling, M. Keim, G. Reuscher, W. Ossau, G. Schmidt, A. Waag, and L. W. Molenkamp, Nature (London) 402, 787 
(1999); Y. Ohno, D. K. Young, B. Beschoten, F.Matsukura, H. Ohno, and D. D. Awschalom, Nature (London) 402, 790 
(1999). 

E. R. Mucciolo, C. Chamon, and C. Marcus, Phys. Rev. Lett. 89, 146802 (2002); Susan K. Watson, R. M. Potok, C. M. 
Marcus, and V. Umansky, Phys. Rev. Lett. 91, 258301 (2003); M. G. Vavilov et ai, Phys. Rev. B 71, 241309(R) (2005). 
5 T. Aono, Phys. Rev. B 67, 155303 (2003) 

" E. Cota, R. Aguado, C.E. Creffield and G. Platero, Nanotechnology 14, 152-156 (2003) 

R. Hanson, L.M.K. Vandersypen, L.H. Willems van Beveren, J.M. Elzerman, LT. Vink and L.P. Kouwenhoven, Phys. Rev. 

B 70, 241304 (2004). 
® P. Recher, E. V. Sukhorukov and D. Loss, Phys. Rev. Lett. 85, 1962 (2000). 
^ M. Blaauboer and C.M.L. Fricot, Phys. Rev. B 71, 041303(R) (2005). 
" G. Platero and R. Aguado, Physics Reports, 395, 1 (2004) and references therein. 

^ W. G. van der Wiel et ai, 'Photon assisted tunneling in quantum dots', in: I.V. Lerner, et al. (Eds.), Strongly Correlated 

Fermions and Bosons in Low-dimensional Disordered Systems, Kluwer Academic Publishers, pp. 43-68 (2002). 
^ Qing-feng Sun, Hong Guo, and Jian Wang, Phys. Rev. Lett. 90, 258301 (2003) 
^ R. Lopez, R. Aguado, G. Platero and C. Tejedor, Phys. Rev. B 64, 075319 (2001) 



10 



J.M. Elzcrman, S. de Pranceschi, D. Goldhaber-Gordon, W.G. van der Wiel and L.P. Kouwenhoven, J. Low Temp. Phys. 
118, 375 (2000) 

15 C. E. Creffield and G. Platero, Phys. Rev. B 65, 113304 (2002); 69, 165312 (2004) 

J.M. Elzerman, R. Hanson, L.H. Willems van Beveren, B. Witkamp, L.M.K. Vandersypen and L.P. Kouwenhoven, Nature 
430, 431 (2004) 

R. Hanson,L.H. Willems van Beveren, I.T. Vink, J.M. Elzerman, W.J.M. Naber, F.H.L. Koppens, L.P. Kouwenhoven and 
L.M.K. Vandersypen, Phys. Rev. Lett. 94, 196802 (2005) 

Ernesto Cota, Ramon Aguado and Gloria Platero, Phys. Rev. Lett. 94, 107202 (2005) 
1" C. A. Stafford and N. S. Wingreen, Phys. Rev. Lett. 76, 1916 (1996). 

B. L. Hazelzet, M. R. Wegewijs, T. H. Stoof, and Yu. V. Nazarov, Phys. Rev. B 63, 165313 (2001). 
T. H. Oosterkamp et al., Nature (London) 395, 873 (1998). 
J. R. Petta et al., Phys. Rev. Lett. 93, 186802 (2004). 
K. Ono and S. Tarucha, Phys Rev. Lett. 92, 256803 (2004) 
F.H.L. Koppens, J. A. Folk, J.M. Elzerman, R. Hanson, L.H. Willems van Beveren, I.T. Vink, H.P. Tranitz, W. Wegscheider, 
L.P. Kouwenhoven and L.M.K. Vandersypen, Science 309, 1346 (2005) 
^5 J.R. Petta et al, Science 309, 2180 (2005) 

K. Blum, Density Matrix Theory and Applications (Plenum Press, NY, 1996) 

Toshimasa Fujisawa, Yasuhiro Tokura and Yoshiro Hirayama, Phys. Rev. B 63, 081304(R) (2001) 

A detailed study of spin relaxation and decoherence in a GaAs quantum dot due to spin-orbit interaction can be found in, 

Vitaly N. Golovach, Alexander Khaetskii, and Daniel Loss, Phys. Rev. Lett. 93, 016601 (2004). 

Rafael Sanchez, Ernesto Cota, Ramon Aguado and Gloria Platero, Physica E (Amsterdam), doi:10.1016/j.physe.2006.03.154 
(2006). 

^° See also, Hans-Andreas Engel and Daniel Loss, Phys. Rev. Lett., 86, 4648 (2001). 



16 



18 



